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SIMULATION  OF  LASER  BEAM  NONUNIFORMITY  EFFECTS 


I.  Introduction: 

In  laser  fusion,  spherical  layered  pellets  are  imploded  through  laser 
induced  ablation  of  the  outer  surface.  The  laser  energy  is  absorbed  near 
the  critical  density  and  the  heat  is  conducted  inward  to  the  cold  ablator 
material.  The  heated  material  expands  outward  and  the  remaining  cold 
material  is  accelerated  inward  as  a  result  of  momentum  conservation.  In 
order  to  achieve  the  densities  and  temperatures  required  for  thermonuclear 
fusion,  this  inward  acceleration  must  be  very  symmetric  over  the  surface 
of  the  pellet.  For  laser  fusion  to  be  practical,  the  laser  energy  must  be 
efficiently  coupled  to  the  pellet.  Implosion  symmetry  and  energy  efficiency 
may  be  conflicting  requirements. 

1  2 

Experiments  and  theory  at  NRL  and  elsewhere  indicate  that  the  most 

efficient  laser  energy  coupling  occurs  at  relatively  low  irradiances 
13  14 

(10  -  10  Watts/  2) .  These  low  irradiances  require  thin-shelled 

cm 

pellets3  (of  thickness  AR)  to  be  accelerated  over  large  distances  (R)  such 

that  R/AR  >  10.  This  places  severe  requirements  on  ablation  pressure 

4 

uniformity  (Ap/p  <  1%)  if  high  gain  is  to  be  achieved.  One  of  the  major 
causes  of  ablation  pressure  nonuniformity  is  nonuniformity  of  the  laser 
beam  intensity.  The  problem  we  address  in  this  report  is  the  relation 
between  nonuniformities  in  the  laser  beam  intensity  and  nonuniformities 
in  the  ablation  pressure  of  thin  planar  targets  as  a  function  of  laser 
intensity  and  scalelength  of  the  inhomogeneity.  These  thin  planar  targets 
model  one  section  of  a  large  radius  thin  shell  pellet. 

An  illustration  of  the  laser-target  interaction  is  given  in  Figure  1. 
The  distance  from  the  ablation  surface,  the  region  of  steepest  density 
Manuacript  submitted  February  19,  1981. 
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gradient  and  where  the  pressure  is  "applied",  to  the  critical  surface, 
where  the  energy  is  absored,  is  D  .  The  scale  length  of  the  inhomogeneity 

dC 

in  laser  intensity  is  \^.  If  the  distance  from  the  ablation  surface 

to  the  critical  surface  is  much  larger  than  the  scalelength  of  this 

inhomogeneity  (i.e.,  D  »X^)  then  lateral  energy  flow  should  smooth  out 

the  inhomogeneity  before  it  reaches  the  ablation  surface. ^  If  D  «  then 

ac  ■* 

the  inhomogeneity  should  imprint  itself  on  the  target  and  the  variation  in 

ablation  pressure  should  be  of  the  same  order  as  the  variation  in  laser 

intensity.  This  inhomogeneity  in  ablation  pressure  will,  in  turn,  affect 

the  target  velocity  and  acceleration  profiles. 

6  7 

Theory  and  experiment  indicate  that  the  distance  from  the  ablation 
surface  to  the  critical  surface  is  directly  proportional  to  the  intensity 
of  the  laser.  Here  we  determine  the  extent  to  which  lateral  energy  flow  can 
smooth  out  laser  inhomogeneities  as  a  function  of  laser  intensity. 

In  the  next  section  we  discuss  our  model.  Section  III  contains  the 
results  cZ  the  simulation  with  the  summary  and  conclusions  in  Section  IV. 
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II.  The  Models 


We  model  the  effect  of  an  inhomogeneous  laser  beam  on  a  thin  (15pm 

8  9 

thick)  plastic  (CH)  foil  with  the  FAST2D  ’  laser-shell  simulation  code. 

This  is  a  fully  two-dimensional  Cartesian  code  with  a  sliding  Eulerian 

grid.  Tt  solves  the  ideal  hydrodynamic  equations  using  the  flux-corrected 

transport  (FCT) ^  algorithms  with  a  two-dimensional  classical  plasma 

thermal  conduction  routine.  The  thermal  conduction  coefficient  has  the 
5/2 

classical  KQT  temperature  dependence  where  KQ  is  the  standard  Braginskii 

11  8  7 /2 

coefficient,  3.06x10  erg/ (cm-sec-ev  '  ). 

The  grid  is  initially  uniform  in  the  region  of  the  foil.  Inside 

(opposite  the  laser)  and  outside  (the  underdense  plasma)  of  this  region 

the  grid  is  stretched.  The  inside  region  is  filled  with  a  very  low 

density  plasma  having  the  same  adiabat  as  the  target  plasma.  On  the 

outside,  beyond  the  critical  surface,  both  the  density  and  pressure  drop 

off  exponentially.  The  system  is  periodic  in  the  transverse  direction. 

All  the  laser  energy  is  assumed  to  be  absorbed  through  resonant 

absorption.  The  laser  flux  is  introduced  as  an  energy  source  in  the 

heat  equation  and  is  deposited  into  the  two  cells  about  the  critical 

density.  The  plasma  is  assumed  to  be  fully  ionized. 

The  initial  density,  pressure  and  temperature  profiles  are  generated 

from  a  one-dimensional,  analytic,  quasi-static  equilibrium  solution  of  a 

12 

set  of  hydrodynamic  equations  that  model  a  laser-plasma  interaction. 

The  solutions  have  been  shown  to  have  provided  an  adequate  steady  state 
with  a  one-dimensional  laser-foil  computer  code. * 3  A  representative 
initial  profile  is  illustrated  in  Figure  2.  The  initial  conditions  are 
assumed  uniform  in  the  transverse  direction  of  the  foil.  A  two-dimensional 
perspective  plot  of  the  density  is  given  in  Figure  3. 


3 


Laser  beam  profiles  are  typically  Gaussian  with  inhomogeneities 
manifesting  themselves  by  a  reduced  intensity  in  the  center  of  the  beam. 
See  Figure  4.  Since  we  are  only  interested  in  the  effect  of  laser  beam 
asymmetries,  we  invoke  the  periodicity  of  the  code  and  model  only  the 
peak-to-peak  variation  in  intensity.  is  the  scale  length  of  the 
inhomogeneity.  The  results  of  a  parameter  variation  in  laser  intensity 
and  scalelength  of  the  inhomogeneity  are  presented  in  the  next  section. 


III.  Results: 

12  13 

Three  different  average  laser  intensities,  <I>  =  5x10  ,  1x10 

13  2 

and  2x10  W/cm  are  used  with  the  scalelength  of  the  inhomogeneity  ranging 

between  100  microns  and  600  microns.  For  all  the  results  presented  here 

the  variation  in  intensity  is  80%  (— §•  =  0.80  .  AI  =  I  -  I  .  )  and 

<I>  max  min 

the  laser  wavelength  is  1.06  pm.  Data  are  shown  at  3  and  6  nanoseconds. 

Figure  5  shows  the  variation  in  the  ablation  pressure  as  a  function  of 

the  scalelength  of  the  inhomogeneity  for  the  three  intensities  after  3 

nanoseconds.  Note  that  the  large  scalelength  variation  (A^  =  600  pm) 

severely  impacts  the  target  for  all  three  intensities.  For  Case  A, 

12  2 

<I>  =  5.0x10  W/cm  ,  lateral  energy  flow  has  no  significant  impact 

for  scalelength  greater  than  100  pm  whereas  for  Case  C,  <I> =  2.0xl013 
2 

W/cm  ,  lateral  flow  gives  rise  to  only  a  5%  pressure  variation  for 
A^  =  500  pm. 

In  Figures  6-9,  we  present  a  detailed  look  at  the  foils  for  Case  B, 

13  2 

<I>  =  1.0x10  W/cm  ,  with  A^  =  100  pm  and  600  pm.  A  perspective  plot 

-9 

of  the  density  at  T  =  6x10  s  for  A^  =  600  pm  is  shown  in  Figure  6.  As 
can  be  seen,  the  ablation  rate  variation  follows  the  laser  intensity 
variation  and  the  foil  becomes  severely  distorted  and  is  beginning  to 
buckle.  Figure  7  illustrates  the  pressure  contours  at  the  same  time. 

Here  the  laser  light  is  coming  from  the  right.  The  lines  are  contours  of 
constant  pressure  in  0.10  increments  counting  from  the  outside  inward. 

The  ablation  surface  and  critical  surface  are  also  indicated.  Note  the 
pockets  of  very  high  pressure  at  the  "edges"  of  the  ablation  surface. 

These  are,  of  course,  the  regions  of  highest  laser  intensity. 

The  density  plot  for  A^  =  100  pm  is  given  in  Figure  8.  Here  there  is 
negligeable  distortion  of  the  foil.  The  corresponding  pressure  contours 


ps* 


are  presented  in  Fig.  9.  In  this  case  the  ablation  surface  is  essentially 
equivalent  to  a  contour  surface.  Lateral  energy  flow  has  provided 
enough  smoothing  to  produce  only  a  0.5%  variation  in  the  ablation  pressure. 

XV.  Summary  and  Conclusion; 

The  numerical  results  are  summarized  in  Figure  10,  where  the  variation 

in  ablation  pressure  is  plotted  as  a  function  of  the  average  laser 

intensity  for  all  six  scalelengths.  This  is  equivalent  to  plotting  the 

variation  in  ablation  pressure  as  a  function  of  the  average  distance 

between  the  ablation  surface  and  the  critical  surface  as  indicated  on  the 

graph.  The  direct  correlation  between  Dac  and  the  laser  intensity  has 

7  6 

also  been  noted  experimentally  and  theoretically. 

Lateral  energy  flow  in  the  form  of  transverse  thermal  conduction 

does  tend  to  smooth  the  effects  of  laser  beam  nonuniformities.  The  degree 

of  smoothing  is  directly  related  to  the  distance  from  the  ablation  surface 

14 

to  the  critical  surface.  To  smooth  non-uniformities  to  the  extent  that 

the  ablation  pressure  variation  is  less  than  1%  requires  that  D,„ ~0  (X  )  . 

This  result  is  somewhat  intuitive  in  that  if  the  scalelengths  of  temperature 

variations  are  comparable,  a  thermal  wave  propagating  transverse  to  the 

laser  beam  should  smooth  out  moderate  temperature  variations  at  about  the 

time  that  the  laser  energy  has  been  conducted  to  the  ablation  surface. 

6 

Detailed  one-dimensional  numerical  results  indicate  that  the  scaling 
for  the  distance  from  the  ablation  surface  to  the  critical  surface  has 
the  following  relation 

D  a  (I/u3) 

3C 

where  w  is  the  frequency  of  the  laser  light.  The  results  presented  here  give 
the  same  scaling  with  intensity  and  indicate  that  laser  intensity  variations  can 
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be  smoothed  so  that  they  have  a  negligeable  impact  on  the  target  as  long 
as  the  scalelength  of  the  inhomogeneity  is  comparable  to  the  distance 
from  the  ablation  surface  to  the  critical  surface. 

These  results  indicate  that  laser  asymmetries  must  be  considered 
in  the  design  of  a  reactor-size  laser  fusion  system.  Employing  state-of- 
art  techniques,  it  is  unlikely  that  the  illumination  uniformity  of  the  target 

plane  would  be  much  better  than  40%.  Using  a  many  beam  system  with  severe 

16 

overlap,  each  beam  would  have  a  spotsize  on  the  order  of  the  target  radius 

(2-3mm) ,  giving  rise  to  intensity  variations  of  the  same  order.  To  smooth 

intensity  variations  with  a  2-3mm  scalelength  would  require  laser  intensities 

on  the  order  of  2-3x10  w/cm  ,  assuming  a  1.06pm  laser.  (Since  the 

density  falls  off  much  more  rapidly  in  spherical  geometry  than  in  slab 

geometry,  the  required  intensity  to  give  D  3mm  is  likely  to  be  even 

ac 

14  2 

greater  than  3x10  w/cm  ).  The  benefits  of  increased  smoothing  must 

then  be  balanced  against  poor  absorption  and  plasma  instabilities  that 

occur  at  high  irradiances.  It  would  seem  that  the  most  likely  means 

to  increase  the  distance  from  the  ablation  surface  to  the  critical 

14,  15 

surface  is  to  reduce  the  frequency  of  the  laser  light,  contrary  to  the 

17 

current  lines  of  research. 

The  results  obtained  here  are  qualitatively  similar  to  the 

experimental  results  obtained  at  NRL  where  opaque  strips  were  used 

18 


to  perturb  the  laser  intensity. 
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interaction,  DflC  is  the  distance  from  the  ablation  surface 
the  scalelength  of  the  laser  intensity  variation. 
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Fig.  2  -  A  semi-log  plot  of  the  initial  pressure,  density 
and  temperature  profiles  of  a  thin  plastic  foil. 
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Fie.  3  -  A  two-dimensional  perspective  plot  of  the  initial  11,888 
density.  The  laser  beam  is  coming  in  from  the  front  with  an 
inhomogenity  scale-length  of  transverse  to  the  direction  of 
propagation. 
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Fig.  5  -  The  variation  in  ablation  pressure  (dp/<p>)  as  a  function  of 
the  scale-length  of  the  inhomogenelty  ('■  )  for  three  different  average 
laser  intensities  (<I>«  5.0  x  10,  1.0 x  1013  and  2.0  x  1013  W/cm"). 
:ll/<I  ^  *0.80  for  ail  three  intensities. 
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Ft*.  7  -  Pressure  contours  for  <I>  =  l.OxlO13  W/ctn  and  \  -  100  after  6  nanoseconds.  The  lines 
are  contours  of  constant  pressure  in  0.10  increments  of  the  maximum  pressure.  The  critical  and  abla 
tion  surfaces  are  also  indicated. 


perspective  density  plot  for  <i>  =  1.0xl01J  W/cm z  and  =  100  nn  after  6  nanoseconds. 
At  this  short  scalelength,  very  little  distortion  is  observed. 
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Pressure  contours  for  <!>  =  1.0  x  10  W/cm  and  X.  =  100  nm  after  6  nanoseconds 


Fig.  10  -  The  variation  in  ablation  pressure  (Ap/<p>)  for  all  6  scalelengths 
plotted  as  a  function  of  the  average  laser  intensity  (<I>)  and  as  a  function 
of  the  average  distance  from  the  ablation  surface  to  the  critical  surface 

(<D  >) 


